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PREFACE 


This  report  summarizes  results  of  a  one  year  study  dealing  with  the  development 
of  analytical  methods  for  predicting  response  of  energy-absorbing  foundations  on  ships 
undergoing  nuclear  weapons  attack.  In  the  first  phase  of  the  study,  a  survey  was  made  of 
the  literature  dealing  with  various  energy-absorbing  configurations  and  the  associated 
dynamic  material  behavior.  n  the  second  phase,  study  was  begun  of  configurations  that 
showed  promise  for  use  in  ships,  namely  side- loaded  tubes,  axially  loaded  tubes,  honeycomb, 
and  beams.  Tests  were  conducted  in  conventional  test  machines  at  various  rates  and  in  a 
drop-table  shock  machine.  Analytical  methods  were  developed  for  side- loaded  tubes  and 
beams. 

This  work  was  sponsored  by  the  Defense  Nuclear  Agency  under  RDT&E  RMSS  Code 
B3440  78464  Y99QAXSF503-11  H2590D. 

The  benefit  of  technical  discussions  with  Lt.  R.  Wade,  Lt.  R.  Elsbemd,  and  Dr. 
Nicholas  Perrone  is  gratefully  acknowledged.  The  computer  analysis  of  the  simulated 
machine- foundation  configuration  under  shock  machine  loading  was  performed  by  Mr.  Koorosh 
Kasraie,  graduate  student  in  the  Department  or  Engineering  Science  and  Mechanics  at  The 
Pennsylvania  State  University. 
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Table  1.  Conversion  factors  for  U.  S.  customary  to  metric 
(SI)  units  of  measurement. 


To  Convert  From 

To 

Multiply  by 

foot 

meter  (m) 

3.048000 

X  E 

-1 

foot-pound- force 

joule  (J) 

1.355818 

inch 

meter  (m) 

2.540000 

X  E 

-2 

kip  (1000  lbf) 

newton  (N) 

4.448222 

X  E 

+  3 

kip/inch2  (ksi) 

kilo  pascal  (kPa) 

6.894757 

X  E 

+  3 

pound- force 

newton  (N) 

4.448222 

pound- force  inch 

newton-meter  (N»m) 

1.129848 

X  E 

-1 

pound- force/ inch 

newton/meter  (N/m) 

1.751268 

X  E 

+  2 

pound- force/ foot 2 

kilo  pascal  (kPa) 

4.788026 

X  E 

-2 

pound-force/inch2  (psi) 

kilo  pascal  (kPa) 

6.894757 

pound-mass 

kilogram  (kg) 

4.535924 

X  E 

-1 

pound-mass-foot2 

(moment  of  inertia) 

kilogram-meter2  (kg-m2) 

4.214011 

X  E 

-2 

pound-mass/foot 3 

kilogram/meter3  (kg/m3) 

1.601846 

X  E 

+  1 

slug 

kilogram  (kg) 

1.459390 

X  E 

+  1 
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H3CEDXMQ  PAOI  BUMkNOT  FILMED 


I .  BACKGROUND 


1.  GENERAL 

This  report  summarizes  results  of  the  initial  phases  of  an  investigation  of  the 
application  of  energy-absorbing  foundations  for  protection  of  equipment  in  ships  under 
nuclear  weapon  attack.  The  results  apply  in  general  to  structures  submitted  to  base 
motion . 

There  are  three  parts  of  the  report.  First  a  summary  is  given  of  the  pertinent 
publications.  Second,  an  outline  is  given  of  analytical  procedures  used  during  the  first 
year  of  the  study.  Third,  experimental  and  analytical  results  are  given  from  the  first 
phase  of  the  present  study  on  side-loaded  tubing,  end-loaded  tubing,  honeycomb,  and  beams. 

With  regard  to  energy- absorbing  foundations,  the  parameters  of  particular  interest 
are  associated  with  material  behavior  and  configuration  characteristics. 


2.  MATERIAL  BEHAVIOR 

The  important  material  properties  are  elasticity,  yield  stress,  ductility,  rate 
sensitivity,  post-yield  stiffness,  and  ultimate  strength. 

The  metals  of  interest  are  steel  and  aluminum.  The  strain  rates  are  associated 
with  gross  structural  motions  of  equipment  on  their  foundations  and  range  up  to  about 
30/sec.  The  elastic  modulus  is  practically  insensitive  to  such  loading  rates.  A  pri¬ 
mary  effect  of  rate  is  the  increase  in  yield  strength.  Hecker  [2]  recently  summarized 
experimental  studies  of  yield  phenomena  in  biaxially  loaded  metals,  giving  277  refer¬ 
ences.  Lindholm  [3]  reviewed  methods  for  conducting  tests  over  the  range  of  strain 

rates.  Steel  is  typically  very  strain  rate  sensitive,  while  aluminum  is  much  less  so. 

•  • 

The  dependence  of  the  dynamic  lower  yield  stress  a(e)  on  the  plastic  strain  rate  c,  for 
strain  rates  up  to  1000/sec  was  represented  by  Cowper  and  Svmonds  [4]  as  the  following 
empirical  expression,  1_ 


o(e)  _ 


1  + 


(1) 


where  a  is  the  static  vield  stress,  D  and  p  arc  material  constants  and  p  mav  be  taken 
o  '  \ 

as  an  integer.  The  values  D  =  40  sec  and  p  =  5  were  deduced  from  data  from  tension 
impact  tests  of  steel  of  Manioine  [5],  For  aluminum  [6]  the  values  are  D  =  6500/sec  and 
p  =  4.  Belsheim  [71  discussed  dynamic  yield  stress  as  related  to  strain  rate  and  Vigness 
[8]  summarized  data  relating  dynamic  yield  to  "de layed-vield  time".  The  time 

delay  between  applied  dynamic  stress  and  resulting  strain  is  another  important  effect 
which  is  accounted  for  to  some  extent  in  the  constitutive  equation  proposed  by  Malvern  [9] 


€ 


1 


a 

—  + 

E 


K(a-ast)n 


(2) 


Here  e  and  a  are  functions  of  time,  E  is  Young's  modulus,  cr  is  the  static  stress  for 
the  same  strain,  and  K  and  n  are  material  parameters.  Cristescu  [10]  suggested  a  consti¬ 
tutive  equation  of  the  type: 

e  =  2.  +  $(o,e)  a  +  f(o,£)  (3) 

where 


V(o,e)  = 


K(e) 


[a-f(e)]  if  a>f(c) 


if  a<f(e) 


(4) 


The  function  $  is  the  measure  of  the  instantaneous  plastic  response.  The  T(o,e)  is 
similar  to  the  second  term  in  equation  (1),  but  it  depends  on  the  overstress  [a-f(E)], 
where  f(e)  is  the  relaxation  boundary,  below  which  the  response  is  elastic  (see  Fig.  1). 

The  relaxation  boundary  and  the  static  stress-strain  curve  may  not  be  the  same.  On  the 
dynamic  stress-strain  curve  in  Fig.  1,  the  maximum  stress  occurs  at  B  and  the  maximum 
strain  at  C.  The  time  delay  between  occurrence  of  maximum  stress  and  maximum  strain  is 
the  relaxation  time. 

The  type  of  law  in  equation  (1)  applies  best  to  perfectly  plastic  materials 
since  it  relates  only  to  increased  yield  stress.  Perrone  [11]  suggested  a  flow  law  in 
Equation  (5]  which  accounts  for  the  strain  hardening  aspects  as  well  as  increased  yield: 

I 

[1  +  f1  +  ^ 

o  (5) 


c  =  f(e) 


(6] 


If  the  strain  hardening  is  a  function  of  strain  rate,  then  equation  (6)  could  be  used 
with  (5). 

Plass  (12]  developed  a  constitutive  equation  for  beam  bending  by  assuming 
plane  sections  remain  plane  after  bending  and  integrating  (2)  over  the  beam  cross-section, 
with  the  result 


k  = 


l 

R(M-Mst)n 


(7) 


Here  k  is  curvature  rate,  M  is  dynamic  moment,  El  is  bending  stiffness,  and  M  is  moment 
obtained  from  the  static  moment -curvature  curve. 
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3. 


CONFIGURATION'  CHARACTERISTICS 

The  important  configuration  characteristics  are  energy-absorbing  efficiency, 
elastic  stiffness,  ultimate  deformation,  directional  properties  and  transmitted  acceler¬ 
ations  . 

3.1  Side-loaded  tubes 

The  large  plastic  response  of  side-loaded  tubes  under  static  loading  has  been 
studied  analytically  and  experimentally  [13,  14].  The  deformation  £  is  related  to  load 
P  and  tube  diameter  d  by 

P  = 

The  solution  ceases  to  be  valid  for  &  >  .  The  load  P„  is  the  load  at  which  the  tube 

Jv  0 

begins  to  deform  as  a  mechanism  after  formation  of  plastic  hinges  at  the  quarter-points. 

The  load  PQ  is  related  to  the  fully  plastic  moment  M0 ,  the  mean  radius  r,  the  tube  thick¬ 
ness  t,  and  the  tube  width  W  by 

4M  a  t2W 
p  -  — fi.  =  J2 _ 

«r  r  (9) 

For  the  dynamic  behavior  of  material  under  uniaxial,  impulsive  loading  Perrone 
[15]  suggested  a  simplification  in  the  analysis,  namely  that  the  initial  strain  rate 
could  be  assumed  to  be  constant  throughout  the  entire  flow  process.  For  a  pulse  loaded 
structure,  it  is  necessary  to  estimate  the  peak  strain  rate  and  assume  the  associated 
stress  to  be  constant. 

For  an  application  of  side-loaded  rings  to  highway  impact  data,  Perrone  [16] 
presented  stress  versus  strain  rate  data  for  10  2<  e  <  400/sec  and  static  load  versus 
deflection  curves  for  tubes  of  three  different  types  of  steel.  These  tubes  were  18 
inches  in  diameter. 

With  regard  to  energy  absorption,  Perrone  assumed  that  the  1"-,,!- deformation 
curve  rose  linearly  from  PQ  at  zero  deflection  to  2P0  at  deflection  equal  to  ir,  so 
that  the  energy  U  absorbed  is  approximately 

l)~3P0r  =  12M0  =  3ootzK  (10) 

The  result  is  surprising,  since  the  energy  absorbed  does  not  depend  on 
radius  of  the  tube.  Rate  of  loading  effects  may  be  taken  into  account  by  increasing  the 
yield  stress  according  to  equation  (1). 

It  should  be  noted  that  in  reference  [14]  it  was  assumed  that  the  tube  was 
being  crushed  between  parallel  plates  so  that  the  contact  point  would  move  out  as  the 
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characteristic  peanut  shape  develops  for  the  tube.  In  reference  [16],  the  loading  appears 
to  be  more  centrally  applied,  which  results  in  a  softer  load-deformation  curve  for  large 
deformations. 
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I  The  maximum  load  carried  should  occur  at  maximum  deformation,  unless  the  tube 

fails  before  that,  failure  could  occur  at  the  plastic  hinges  for  brittle  metals.  The 

f  suitability  of  the  material  can  best  be  determined  by  tests  of  the  tubing  itself. 

1. 

:  3.2  Monevco'ii 

!  Hone. Co"  >-  received  considerable  attention  as  an  energy-absorbing,  cushion¬ 

ing  material  fr  ;  'e.tian  during  air  drops  of  equipment,  moon  landing  of  space  craft 
and  for  many  other  app 1 i cat i ons .  An  excellent  tutorial  paper  on  cushioning  for  aerial 

•  delivery  was  given  by  IT  pson  and  Ripperger  [17]. 

The  materials  which  have  been  tested  most  are  paper  and  aluminum  honeycombs. 

For  air  drops,  the  ideal  load-deformation  curve  is  flat.  One  measure  of  the  efficiency 
for  energy  absorption  is 

Efficiency  =  xl00%  _ 

m  *• 


where  a  =  maximum  stress  and  d  is  maximum  deformation.  Thus  a  rectangular  stress-defor- 
m 

mation  curve  would  represent  the  ideal  efficiency  of  100%.  Paper  honeycomb  has  a  nearly 
ideal  stress-strain  curve  up  to  about  70%  strain,  above  which  it  stiffens  as  it  begins 
to  bottom  out. 

Aluminum  honeycomb  stress-strain  curves  typically  show  an  initial  sharp  stress 
peak,  called  the  compressive  strength,  at  which  buckling  begins.  As  the  strain  increases 
further,  the  stress  drops  to  a  nearly  constant  level,  called  the  crush  strength,  which 
is  usually  about  65  to  70  per  cent  of  the  compressive  strength.  Data  sheets  are  available 
which  list  compressive  strength,  crush  strength,  shear  strength  in  two  directions,  and 
beam  shear  modulus.  A  typical  aluminum  honeycomb  may  be  specified  as  3.7  -  3/8  -  .0025 
which  means  successively:  density  in  pounds/ft3,  cell  size  in  inches,  and  metal  gage 
in  inches. 


With  regard  to  analysis,  McFarland  [18]  discussed  two  possible  modes  of  failure 


where 


t  =  cell  wall  thickness 
S  =  cell  minor  diameter 
D  =  3/3^  =  width  of  cell  wall 

=  width  of  basic  panel  element  [D/4<Pw<D/2] 

The  deformation  occurs  in  an  accordian-like  pattern,  with  Pw  the  length  of  one  pleat, 
measured  along  the  slope,  from  peak  to  valley.  For  shear  deformation,  using  q^  as  yield 
stress  for  shear. 


U  =  0.433  q  tDP 
s  "o  w 

1 . 155q  t 

F  =  - -£_ 

3  s 

The  total  energy  U  and  total  force  F  are 

U  =  U  +  u 
c  s 


(14) 

(15) 


(16) 


F  =  F  +  F 
c  s 


(17) 


In  addition  it  is  noted  in  [18]  that  if  the  ratio  t/S  >  0.004,  the  gross  shear 
mode  failure  will  occur,  rather  than  the  crushing  mode.  This  is  undesirable  since  the 
energy  absorption  in  the  shear  mode  is  less  than  in  the  crushing  mode  for  the  same  t/S 
ratio. 

For  aluminum  honeycomb,  the  initial  sharp  peak  in  the  stress-strain  curve  is 
undesirable  and  may  be  eliminated  [19]  by  a  slight  precrushing  or  pre-dimpling  of  the 
surfaces . 

3.3  Axially  Loaded  Crushed  Tubes 

Axially  loaded  tubing  has  been  investigated  in  a  crushing  mode  [20,  21],  a 
"frangible"  mode [22] ,  and  an  inverting  mode  [23], 

The  axially  loaded  tubing  in  a  crushing  mode  deforms  in  an  accordian-like 
fashion.  The  load-deformation  curve  has  a  large  initial  peak  and  then  oscillates  about 
a  fairly  flat  average.  The  oscillations  are  due  to  the  successive  formation  and  bottom¬ 
ing  out  of  new  pleats. 

Shaw  [24]  studied  this  type  of  deformation  and  attempted  to  predict  the  length 
h  of  tube  involved  in  one  pleat.  An  Euler  instability  is  assumed  to  be  involved.  The 
buckling  load  for  a  slender  column  of  height  h  with  hinged  ends  constrained  to  remain 
vertically  aligned  is 
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p  -  i!ei 

h2 


If,  further,  the  instability  occurs  when  the  tube  just  becomes  plastic,  then 


(17) 


P  =  TTdtO  (18) 

Then,  using  I  =  ~-frdt3,  eliminating  P  gives 

h  _  /  tt  5  E~  "  1 

1  12°o  (19) 

Shaw  said  this  agreed  well  with  the  observed  pattern  with  a  mandrel  in  place,  which  aided 
in  stabilizing  the  formation  of  the  pleats.  The  load  was  related  to  deformation  x  by 


P  =  lrdtao  (5f  +  2 \/2\  ■  1  1  j  (2 

The  minimum  value  of  x  for  which  this  equation  holds  is  xq,  the  displacement  at  which  the 
plastic  hinge  first  develops.  The  value  of  xq  is 


1  + 


The  approximate  outside  diameter  of  the  flange  formed  in  one  pleat  is 


(21) 


D  =  d  +  h  -  t  (22) 

The  predicted  value  from  these  equations  agreed  well  with  static  test  results,  but  not 
with  dynamic. 

3.4  Frangible  Tubes 

The  frangible  tube  as  an  energy- absorbing  device  was  studied  by  McGchce  [25,26]. 
In  the  fragmenting-tube  process,  an  axial  load  is  applied  to  one  end  of  the  tube  while 
the  other  end  is  pressed  over  a  die.  The  die  is  shaped  so  that  the  portion  of  the  tube 
in  contact  with  the  die  is  split  into  segments  and  the  segments  are  broken  into  small 
fragments.  A  fluctuating  force  is  developed,  but  the  average  force  is  approximately 
constant . 

An  empirical  equation  for  the  fragmenting  stress  in  ksi  for  2024-T5  aluminum- 
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alloy  tubing  is  given  [26]  as 


a 


f 


0.7  -  t/r 


(23) 


where  a ^  =  axial  force/area 

t  =  tube  thickness 
D.  =  inner  tube  diameter 

l 

r  =  forming  radius  of  die 

The  range  of  t/r  studied  was  from  0.333  to  0.644.  For  t/r  less  than  0.25  the  tube  seg¬ 
ments  split  and  roll,  but  do  not  fragment.  It  was  found  for  tubing  of  uniform  thickness 
that  the  initial  force  was  higher  than  that  required  after  fragmenting  began.  To  reduce 
this  peak,  a  14°  taper  was  used  on  the  tubing  thickness  at  the  die  end  of  the  tubing. 

On  page  ,  the  energy- absorbing  capability  of  four  metals  studied  as  frag¬ 
menting  tubes  in  reference  [26]  is  compared  with  that  of  aluminum  honeycomb  [27]  and 
cellular  aluminum  [28],  both  in  a  crushing  mode. 

3.5  Inverting  Tubes 

The  axially  loaded  tubing  which  is  inverted  is  first  flared  on  one  end.  Then 
the  flared  end  is  rigidly  clamped  and  the  tube  is  turned  inside-out  by  pushing  the  tube 
through  the  flared  end.  Kroell  [23]  gave  many  experimental  results  for  3003  aluminum, 
either  annealed  or  half-hard,  and  a  few  for  mild  steel.  The  load-deflection  curves  are 
flat.  The  specific  energy  for  3003-H14  aluminum  varies  from  2000  to  7000  ft-lb/lb  for 
tubing  with  thickness  ranging  from  .015  to  .060  inches  and  mean  diameter  varying  from 
0.8  to  3.0  inches. 

3.6  Energy  Absorbed  by  Cutting  Metal 

In  order  to  increase  the  axial  force,  Kroell  [23]  added  "drag"  elements  around 
the  outside  of  the  tube  which  scratched  the  tubing.  The  total  energy  lost  was  up  to  2.37 
times  that  lost  due  to  inverting  only. 

Shaw  [24]  discussed  an  energy-absorbing  device  in  which  a  circular  machine  tool 
is  designed  to  simultaneously  skin  and  draw  a  round  bar.  For  mild  steel  about  300,000  in. 
lb.  of  energy  is  required  to  produce  one  cubic  inch  of  chips;  for  aluminum  about  150,000 
in. lb.  is  required.  To  save  weight,  the  tube  was  made  thinner  and  then  the  crushing  mode 
occurred  before  the  skinning  could  develop  properly.  Thus,  the  weight  of  the  tubing  used 
for  skinning, for  the  same  application  to  automobile  bumpers,  was  greater  than  that  used  in 
an  alternate  design  in  which  the  tube  was  crushed  axially,  so  it  appears  that  the  tube 
crushed  under  axially  loading  absorbs  more  energy  per  pound  of  material. 

3.7  Beams 

Beams  are  used  to  support  equipment  on  ships,  but  they  arc  usually  designed  to 
behave  elastically;  therefore,  their  energy-absorbing  capacity  is  not  used.  Further,  the 


15 


need  to  design  elastically  often  leads  to  the  choice  of  higher  yield  steels,  which  tend 
to  be  more  brittle.  Elastic  design  must  be  used  where  final  position  and  alignment  must 
be  maintained,  but  there  are  situations  where  permanent  deformations  would  not  interfere 
with  the  mission  of  the  structure.  Therefore,  the  ability  to  perform  inelastic 
analysis  for  ship  shock  needs  to  be  further  developed. 

Much  has  been  published  on  the  dynamic  plastic  behavior  of  structures.  Baker 
[29]  summarized  approximate  techniques.  Jones  [30]  presented  a  literature  review  of  the 
general  topic  and  referred  to  the  following  literature  surveys:  Goldsmith  [31]  on  beams, 
Rawlings  [32]  on  beams  and  frames,  Symonds  [33]  on  beams,  Johnson  [34]  on  beams,  frames 
and  plates,  Jones  et  al.,  [35]  on  plates  and  shells,  and  Lee  [36]  on  beams,  plates  and 
shells.  Krajcinovic  [37]  published  a  survey  of  exact  solutions  for  the  dynamic  rigid- 
plastic  behavior  of  various  perfectly  plastic  structures  and  also  discussed  some  of  the 
bounding  methods. 

While  many  papers  have  been  published,  their  application  is  limited  because  of 
various  assumptions,  such  as:  the  material  is  rigid-perfectly  plastic,  the  load  increases 
monotonically,  the  load  is  impulsive,  the  deflection  curve  assumes  a  certain  shape,  shear 
effects  on  yielding  and/or  on  deformation  may  be  ignored,  plastic  hinges  remain  a  constant 
length,  etc.  A  more  general  approach  may  be  achieved  through  application  of  numerical 
methods  and  finite  element  techniques,  such  as  that  of  Witmer,  Balmer,  Leech,  and  Pi  an 
[38].  They  use  layers  to  represent  the  cross-section  of  the  beam,  plate,  or  shell,  lump 
the  mass  at  discrete  points,  and  do  a  time-wise  integration.  This  tends  to  be  expensive 
on  the  computer,  because  of  the  number  of  nodes  and  layers  required  and  because  the  time 
interval  must  be  small  to  achieve  stability  of  the  numerical  solution.  Improvements  have 
been  discussed  by  Leech,  Witmer,  and  Pian  [39],  Morino,  Leech,  and  Witmer  [40], and  Wu  and 
Witmer  [41,  4_j. 

The  writer  and  associates  have  studied  beams  bent  slightly  into  the  plastic 
range  due  to  dynamic  loading  like  that  occurring  in  ship  shock.  Vogel  [44]  developed  a 
computer  approach  to  analyze  beams  tested  on  the  Barry  shock  machine  at  the  Naval  Research 
Laboratory.  He  accounted  for  large  deflections  and  used  the  special  form  of  constitutive 
equation  for  beam  bending  given  in  liquation  (7).  Brown  [45]  suggested  values  of  the 
parameters  in  the  same  equat ion  based  on  shock  machine  lists  of  steel  beams  having  varying 
amounts  of  initial  cold  work.  Prick  [46]  carried  out  tests  or.  copper-nickel  and  steel 
tubing  in  bending. 

Tests  of  cantilever  beams  carrying  heavy  tip  masses  were  made  on  the  Floating 
Shock  Platform  using  315.7  beams,  with  results  reported  bv  Mutt  et.al. ,  [47].  Later,  a 
three-span  hollow  shaft  of  3.125"  0.  b.  was  a  1 s.  . ested  and  experimental  displacement, 
velocity, and  strain  versus  rime  curves  were  presented  [48].  I'i splucement  versus  time  was 
predicted  using  a  finite-difference  method  and  an  e last i c-perfect 1 y  plastic  moment- 
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curvature  relationship.  Bounds  of  peak  displacement  were  estimated  in  one  case  using  an 
energy  technique  and  an  input  initial  velocity  distribution. 

3.8  Summary,  Configurations 

While  there  has  been  an  effort  to  develop  analysis  procedures  for  elastic-plastic 
beams  under  ship  shock  loading,  there  has  been  little  experience  with  design  of  beams 
deformed  into  the  plastic  range  on  ships.  Other  configurations  have  received  little  atten¬ 
tion  as  energy  absorbing  foundations  or  snubbers  on  ships.  Many  configurations  have  been 
discussed  in  the  guide  for  selection  of  elastic  shock  mounts  prepared  by  Burns  [48],  Typi¬ 
cal  non-linear,  elastic, load-deflection  curves  are  given  for  helical  springs,  pneumatic 
cylinders,  hydraulic  cylinders,  airbags  and  rubber  springs.  Some  shapes  are  presented 
which  might  be  used  as  elastic-plastic  springs,  including  cantilever  beams,  hinged-hinged 
beams,  side- loaded  tubes,  and  U-shaped  mounts.  For  these  four  shapes  equations  are  pre¬ 
sented  for  each  of  the  three  principal  axes  for  the  following:  elastic  spring  constant, 
elastic  deflection,  maximum  moment  and  point  of  maximum  stress,  load  at  which  yield  stress 
is  reached  at  the  most  highly  stressed  point,  and  load  at  which  most  highly  stressed  sec¬ 
tion  becomes  fully  plastic.  Thus  the  guide  is  useful  for  designing  elastic-plastic  shock 
mounts  deformed  slightly  into  the  plastic  range. 

The  use  of  energy-absorbing  configurations  for  collision  protection  of  ships 
was  surveyed  by  Jones  [49].  He  compared  the  energy-absorbing  capability  of  foam-filled 
honeycomb  with  that  of  deck  plating.  A  nest  of  tubes  which  would  be  crushed  axially  was 
also  suggested  and  some  approximate  formulas  were  given  for  energy-absorbing  capabil¬ 
ities  of  the  configurations.  While  the  problem  is  different  from  that  of  the  present 
study,  some  of  the  information  on  configurations  is  applicable. 

When  considering  configurations  for  protection  of  internal  ship  or  submarine 
equipment  against  shock,  the  fact  that  the  shock  may  occur  in  any  direction  must  be  con¬ 
sidered.  Axially  loaded  tubes  as  shock  snubbers  oriented  in  three  principal  directions 
may  be  satisfactory,  provided  the  tubes  are  guided  on  mandrels  to  assure  that  the  tubes 
deform  axially,  but  developing  such  a  system  may  be  difficult.  Side-loaded  tubes  and 
honeycomb  have  an  advantage  because  of  their  directional  properties,  so  these  configur¬ 
ations  have  been  considered  first  .n  the  present  study. 

Energy- absorbing  capability  of  various  configurations  were  compared  by  I'errone 
[50].  His  summary,  along  witli  data  from  reference  [J<>]  is  given  in  Table  J. 
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Table  2.  Energy  Absorbing  Systems 


Specific  Energy 


M  st  em 

ft. lb. 
lb.  X 

103 

Comments 

Inverting  tubes  [5<i] 

2  to  4 

Aluminum,  3"  0.  D.,  .035"  wall 
thickness 

Cellular  aluminum  [50] 

1  to  7 

Aluminum  honeycomb  [50] 

1  to  8 

Strain  dependent,  ranges  up  to  70 
percent  strain 

Foam-filled  honeycomb  [50] 

3  to  8 

Side-loaded  tube  [50] 

0.5  to 

1.5 

Side-loaded  tube  (Sandwich)  [50] 

1.5  to 

5 

Frame  assembly  [50] 

0.3  to 

0.5 

Frangible  tube  [26] 

AZ  31  B  Magnesium  alloy 

14 

2024-Ts  Aluminum  alloy 

30 

7075-T6  Aluminum  alloy 

33.5 

AISI  4130  Steel  tubing 

38.5 

Cold  water  quenched 

Axially  loaded  crushed  tube  [Present 
Study] 

Annealed  1018  steel 
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2"  0.  D.,  .0625"  wall  thickness 


II.  THEORETICAL  APPROACH 


4.  MACHINE-FOUNDATION-SNUBBER  ARRANGEMENT 

The  goal  of  the  analysis  is  to  be  able  to  predict  the  response  of  a  machine- 
foundation-snubber  combination  to  a  base  motion  input.  If  the  base  motion  is  translational, 
the  system  may  be  represented  by  the  one-degree-of-freedom  system  shown  in  Figure  2(a). 
During  normal  operation  the  machine  may  be  supported  on  a  linear  spring  and  dashpot  to 
attenuate  steady  vibration.  Then,  during  a  severe  shock  applied  as  base  motion  at  the 
support,  the  shock  snubber  comes  into  contact  with  the  mass.  A  representation  of  the  load- 

deflection  behavior  of  the  snubber  is  shown  in  Figure  2(b).  The  idealized  curve  used  in 

the  first  phase  of  the  analysis  is  assumed  to  have  an  elastic  portion  of  stiffness  kg  . 
After  a  yield  load  is  reached,  the  slope  of  the  curve  changes  to  k^.  When  unloading 
begins,  it  is  assumed  to  occur  elastically,  so  the  stiffness  is  again  k^.  The  differen¬ 
tial  equations  of  motion  are 

my^  +  k(y^  -  yQ)  +  c(y^  -  yQ)  =  -  L(t)  -  mg  (y^  -  yQ)>  a  (24) 

k(yQ  -  +  c(yQ  -  y t)  =  L(t)  +  R(t)  (y^  -  yQ)>  a  (25) 


For  (yj  -  y  )  <  a,  there  is  no  contact  with  the  snubber  and  L(t)  =  o. 

Equation  (25)  was  used  to  predict  the  response  of  a  mass  supported  on  a  set  of 
side-loaded  tubes.  The  clearance,  a  ,  was  taken  as  zero.  The  yield  load  and  the  stiff¬ 
nesses  k0  and  k^  were  obtained  from  experimental  data.  The  solution  was  developed  using 
a  finite  difference  approach  with  central  differences. 

5.  VISCOPLASTIC  BEAM  RESPONSE. 

The  beam  analyzed  was  rigidly  clamped  at  each  end  such  that  the  supports  could 
neither  rotate  nor  move  axially  (Figure  3).  The  support  was  submitted  to  a  vertical  base 
motion  in  a  shock  machine.  The  beam  carried  a  heavy  central  mass.  The  maximum  deflec¬ 
tions  were  about  five  times  the  beams  thickness.  Both  bending  and  axial  deformation  were 
taken  into  account.  The  procedure  used  was  that  developed  by  Witmer  et  al.,  in  reference 
[38],  The  beam  thickness  is  represented  by  layers,  with  each  layer  having  the  constitu¬ 
tive  relationship  of  equation  (2)  with  n  =  1,  or 


e  =  |  +  K(o  -  Ost) 


(26) 


The  procedure  is  outlined  in  a  doctoral  thesis  by  S.  J.  Yim[51],  so  the  details  are  not  re¬ 
peated  here.  Some  significant  results  of  the  study  are  presented  in  the  Section  IV. 


III.  EXPERIMENTAL  PROCEDURE  AND  APPARATUS 

The  purpose  of  the  experiments  was  two-fold.  First,  load-deformation  curves 
were  found  at  various  loading  rates.  Second,  the  attenuation  produced  by  using  energy¬ 
absorbing  snubbers  was  demonstrated. 

The  procedure  was  to  determine  the  material  and  configuration  behavior  at  low 
and  moderate  loading  rates  in  a  conventional  Tinius-Olsen  testing  machine.  An  apparatus 
was  then  designed  so  that  the  configuration  could  be  tested  in  a  simulated  machine-snub¬ 
ber  arrangement  at  a  high  rate  in  a  drop-table  shock  machine. 

The  instrumentation  used  to  obtain  load-deformation  records  with  the  Tinius- 
Olsen  test  machine  was  that  incorporated  in  the  machine.  The  maximum  deformation  rate 
obtainable  is  20  inches/minute. 

A  sketch  of  the  IMPAC  shock  machine  is  shown  in  Figure  4.  The  machine  consists 
of  a  heavy,  220  pound  table  on  which  the  specimens  are  mounted.  The  table  is  raised  to  a 
specified  height,  from  which  it  is  dropped.  The  impact  is  controlled  by  a  spring-like 
"programmer"  on  the  bottom  of  the  table.  After  the  programmer  contacts  the  rigid  base, 
the  table  undergoes  a  severe  negative  acceleration.  The  acceleration-time  curve  is  approx¬ 
imately  a  half-sine  wave  with  durations  of  5  to  25  milliseconds  obtainable  and  maximum 
accelerations  of  lOOg.  After  rebound,  the  table  moves  freely  upward,  guided  by  the  side 
rails.  At  the  top  of  the  rebound,  air  brakes  are  automatically  activated  which  press 
against  the  rails  and  stop  the  table,  preventing  a  second  impact. 

A  simulated  machine-snubber  configuration  on  the  shock  machine  table  is  shown 
in  Figure  5.  The  machine  is  simulated  by  two  steel  plates,  each  8"  x  1/2"  x  8-3/4".  The 
plates  were  bolted  together  to  form  a  single  mass,  and  the  total  weight  was  20  pounds. 

The  snubbers  shown  are  four  side- loaded  ri^es  in  a  symmetrical  arrangement.  The  plate  and 
rings  were  kept  in  place  by  tape  which  held  the  configuration  in  position  during  the  free 
fall,  but  slackened  and  allowed  the  plate  to  separate  from  the  rings  upon  rebound. 

Three  accelerometers  were  attached  to  the  top  of  the  plate,  at  locations  A,  B 

and  C  in  Figure  5.  Bv  this  arrangement  the  vertical  acceleration  of  the  mass  was  measured, 

as  well  as  rotational  acceleration  about  two  perpendicular  principal  axes  through  the 

center  of  the  top  surface  of  the  mass.  A  fourth  accelerometer  was  glued  to  the  top  of  the 

shock  machine  table,  tv)  measure  input  acceleration  v  .  Two  accelerometers  were  Kistler 

'  o 

808A  with  Kistler  545A  charge  amplifiers.  The  other  two  accelerometers  were  Endevco  2222B 
with  Endevco  2750  charge  amplifiers. 

lor  the  beam  tests,  the  two  Ki:  tier  accelerometers  were  used  to  measure  vertical 
acceleration,  one  on  the  central  mass  and  one  on  the  shock  machine  table.  The  unsupported 
length  of  the  beam  was  8"  (Figure  3).  Strain  gauges  were  located  on  top  and  bottom  of  the 
beam  at  distances  of  1/2",  1”,  2",  and  3"  from  the  support.  The  beam  arrangement  was  also 


tested  statically  in  the  Tinius-Olsen  test  machine.  The  strain  gauges  were  Micromeasure¬ 
ment  EP-08-125BB-120. 
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IV.  RESULTS 


Tests  results  were  obtained  for  side- loaded  tubes,  end- loaded  tubes,  honeycomb, 
and  beams.  Analysis  was  carried  out  for  side- loaded  tubes  and  beams. 

6.  SIDE -LOADED  TUBES 

6.1  Test  Results 

The  tubes  tested  were  all  obtained  from  one  piece  of  round  seamless  mechanical 
tubing  of  cold  drawn  carbon  steel,  with  0.10  -  0.25%  carbon.  The  wall  thickness  was  16 

gauge  or  0.065  inches,  and  the  outside  diameter  was  2  inches.  Some  of  the  tubing  was 

annealed,  after  it  was  found  that  the  as  eceived  tended  to  crack  at  the  plastic  hinges 
at  large  deformations  (see  Figure  6(a)). 

Small  tensile  specimens  were  cut  from  the  tubing  with  the  axis  of  the  specimens 
in  the  direction  of  the  axis  of  the  tubing.  The  specimens  had  a  slight  curvature  and  the 
ends  were  flattened  in  the  grips,  but  this  was  thought  to  have  little  effect  on  the  ten¬ 
sile  properties.  The  results  are  shown  in  Figure  7.  The  annealed  material  had  a  definite 
yield  at  44,000  psi. 

Specimens  2"  long  were  cut  from  the  as-received  tubing  in  a  non-ideal,  but  prac¬ 
tical,  way  by  using  a  plumber's  tube-cutter.  This  produced  a  v-shaped  groove.  The  tubing 
was  tested  in  the  Tinius-Olsen  machine  at  varying  rates.  Plastic  hinges  tend  to  develop, 
first  on  the  top  and  bottom  and  then  on  the  sides  of  the  tubing,  so  finally  there  are  four 
hinges  located  at  90°  intervals  around  the  tube  circumference  (see  Figure  6(b)).  All  the 
as-received  tubes  developed  cracks  at  about  one  inch  deformation.  On  Figure  8,  a  load- 
deformation  curve  is  shown  at  a  loading  rate  of  20"/min.  A  sudden,  temporary,  loss  of 
stiffness  occurred  when  the  cracks  developed.  This  tubing  was  judged  to  be  unacceptable 
in  the  as-received  state  and  all  subsequent  tests  were  done  on  annealed  tubing.  Test 
results  for  annealed  tubing  at  two  rates  are  shown  in  Figure  8.  The  varying  slopes  of  the 
initial  elastic  portion  of  the  curves  were  due  in  part  to  the  raised  edges  produced  by  the 
tube-cutter.  Specimens  for  later  tests  were  cut  and  had  edges  machined  on  a  lathe. 

In  the  shock  machine,  four  1/2"  long  tubes  were  used,  rather  than  one  2"  long 
tube.  Test  results  for  four  side-loaded  annealed  tubes  are  shown  in  Figure  9,  for  two 
different  sets  of  tubes,  to  check  repeatability,  at  a  loading  rate  of  0.05'Vmin.  The  load- 
deformation  curve  for  four  1/2"  long  tubes  was  practically  identical  to  that  of  one  2" 
long  tube,  at  the  same  loading  rate. 

In  Figure  10  the  loaa  versus  deformation  curves  of  four  1/2"  tubes  are  shown  as 
obtained  from  the  linius-Olsen  machine  at  0.0S"/min.  and  0.3"/min.  Also  shown  are  the 
load-deformation  curves  obtained  on  the  shock  machine  due  to  two  successive  15"  drops. 

As  expected,  because  of  rate  effects,  there  is  an  increase  in  yield  load  as  loading  rate 


The  deformations  from  shock  machine  data  were  cbtained  by  double  integration  of 
(y  -  yQ) .  The  dynamic  loads  are  (my  -  mg).  For  the  data  plotted,  the  angular  acceler¬ 
ation  of  the  mass  was  small  enough  to  be  neglected.  In  Figure  11,  the  load-deformation 
curves  for  the  successive  drops  are  plotted  with  the  curve  for  the  second  drop  shifted  to 
the  origin.  Note  that  the  load-deformations  curve  is  nearly  flat  for  the  first  drop,  but 
the  configuration  is  stiffer  during  the  second  drop. 

In  anothe“r  set  of  shock  machine  tests,  a  set  of  four  1/2"  long  annealed  tubes 
was  submitted  to  eight  successive  10"  drops  and  then  a  15"  drop.  Resulting  load-deforma¬ 
tion  curves  for  Drops  1  through  4  and  Drop  9  are  shown  in  Figure  12  and  compared  with  a 
curve  at  0.05"/min.  deformation  rate  on  the  Tinius-Olsen. 

The  shifted  curves  for  Drops  1  through  4  are  shown  in  Figure  13.  Less  energy  is 
absorbed  on  each  successive  drop. 

Transmitted  acceleration  is  compared  with  base  acceleration  in  Figure  14  for  the 
first  four  10"  drops.  For  the  first  drop  the  maximum  transmitted  acceleration,  v  ,  is 
about  26  g  compared  to  a  peak  base  acceleration,  y  ,  of  45  g.  However,  by  the  fourth  drop, 
the  tubes  furnished  no  attenuation. 

In  Figure  15,  the  transmitted  acceleration  is  compared  with  input  acceleration 
for  the  9th  drop.  Here  the  peak  input  acceleration  is  58  g  and  the  maximum  transmitted 
is  67  g. 

The  relative  velocity  versus  deformation  for  the  four  successive  10"  drops  is 
shown  in  Figure  16.  These  appear  to  be  almost  semi-circular  during  the  compressive  phase 
for  the  rings.  The  calculated  rebound  (positive)  velocity  is  shown  and  seems  to  vary  only 
slightly  on  successive  drops. 

Load  rate  versus  deformation  for  the  first  two  successive  10"  drops  are  shown  in 
Figure  17.  The  load  rate  is  practically  zero  except  at  the  beginning  and  the  end  of  the 
deformation  stroke.  This  may  have  some  significance  in  developing  an  approximate  analysis 
procedure. 

6.2  Theoretical  Results 

In  Figure  18,  predicted  load-deformation  curves  are  compared  with  those  from 
measurements.  To  pre  lict  the  results,  an  equation  based  on  (241  was  used,  with  y  = 


my  +  L  (t )  =  -mg  (27) 

or  myr  ♦  L (t )  =  -mg  -  my 

The  Lit)  curve  was  a  bi-linear  curve,  as  in  Figure  2(b),  which  was  estimated 
from  test  data  from  the  linius-Olsen  machine,  but  with  the  yield  increased  to  account  for 
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rate  effects.  The  input  was  the  measured  y  .  Agreement  is  good  except  for  the  unloading 
curve  and  the  residual  deformation  for  Drop  2. 

Predicted  acceleration  versus  time  is  given  in  Figure  19  for  10"  Drop  1  and  in 
Figure  20  for  10"  Drop  2.  In  Figure  20,  the  effect  is  studied  of  doubling  the  theoretical 
elastic  stiffness  of  the  tubing.  The  kg  =  8300  lb. /in.  produces  the  better  fit. 

A  comparison  of  theoretical  and  experimental  load-deformation  curves  for  the 
first  of  two  successive  IS"  drops  is  shown  in  Figure  21.  The  simple  representation  is 
again  fairly  adequate  except  during  unloading,  where  the  predicted  residual  deformation  is 
in  error  by  about  8%. 

In  general,  predictions  are  good  for  side-loaded  tubing,  but  could  be  improved 
by  incorporating  a  constitutive  equation  of  the  form  of  equation  (26) . 

7.  AXIALLY  LOADED  CRUSHED  TUBES 

Several  tests  were  run  at  moderate  loading  rates  in  the  Tinius-Olsen  machine  of 
the  2"  diameter  steel  tubing  under  axial  loading.  A  typical  load-deformation  curve  for 
an  as-received,  2"  long  tube  is  shown  in  Figure  22,  in  comparison  with  that  for  an 
annealed  tube  of  the  same  size.  The  ends  of  the  as-received  tube  were  cut  with  a  tube- 
cutter,  while  those  of  the  annealed  tube  were  machined  square.  The  end  conditions  affect 
the  initial  shape  of  the  curve,  until  the  end  starts  to  curl  in  and  the  first  pleat  begins 
to  form.  The  deformed,  annealed  tube  is  shown  in  Figure  6(c),  showing  two  well-formed 
pleats.  The  pleats  develop  successively  and  the  rise  and  fall  of  the  load-deformation 
curve  is  associated  with  the  formation  of  successive  peaks. 

8.  HONEYCOMB 

Experiments  were  performed  on  aluminum,  5052  alloy  honeycomb  donated  by  Ameri¬ 
can  Cyanamid  Corporation.  The  honeycomb  was  military  grade,  density  3.7  lb. /ft?,  3/8" 
core,  0.0025  metal  gage,  and  5/8"  deep.  The  crush  strength  is  published  as  290  psi  and 
compressive  strength  as  200  psi.  A  quasi-static  load-deformation  curve  from  the  Tinius- 
Olsen  machine  of  a  2"  x  2"  specimen  is  shown  in  Figure  23.  There  were  thirty  complete 
cells.  A  photograph  of  a  typical  deformed  specimen  is  shown  in  Figure  6(d).  The  inner  cells 
deformed  by  the  formation  of  successive  pleats,  much  as  the  axially  crushed  tube  shown  in 
Figure  6(c).  Three  pleats  can  be  seen  in  Figure  6(d).  The  outer  walls  buckle  in  a  mode 
involving  the  total  depth  of  the  specimen.  The  initial  peak  or  crush  strength  is  higher 
than  4  in2  x  290  psi  expected  from  the  advertised  value  but  the  advertised  compressive 
strength  of  4in2  X  200  psi  is  close  to  the  average  after  the  first  peak. 

The  initial  high  peak  in  the  load-deformation  curve  is  undesirable  in  a  shock 
snubber,  since  it  leads  to  a  high  transmitted  acceleration.  To  avoid  this  peak  in  the 
shock  machine  tests,  it  was  decided  to  pre-crush  the  honeycomb  in  the  Tinius-Olsen  machine. 
In  Figure  24,  a  load-deformation  curve  is  shown  during  a  quasi-static  pre-crush  of  a 


2"  x  1-3/4"  honeycomb  specimen  with  twenty-five  complete  cells.  Here  the  pre-crush  was 
terminated  prematurely  after  about  0.015"  deformation,  and  an  undesirable  peak  of  1200  lb. 
occurred  on  re-loading.  The  best  procedure  is  to  continue  the  pre-crush  until  the  load 
drops  to  the  average  level  at  about  650  psi.  The  subsequent  oscillation  in  the  load- 
deformation  curve  is  assumed  to  be  due  to  the  successive  formation  of  pleats.  That  they 
formed  successively  and  not  simultaneously  could  not  be  observed  on  the  present  tests, 
because  the  pleats  formed  only  on  interior  cells. 

After  the  pre-crushing  procedure  was  established  a  2"  x  2"  specimen  with  thirty 
complete  cells  was  pre-crushed  and  then  used  as  the  snubber  in  the  machine-snubber  con¬ 
figuration  in  the  shock  machine.  The  machine  was  simulated  by  the  same  plates,  weighing 
20  lb.,  shown  in  Figure  5.  The  honeycomb  was  substituted  for  the  side-loaded  rings,  but 
the  single  specimen  of  honeycomb  was  positioned  under  the  center  of  the  plates.  In  Figure 
25,  the  load-deformation  curve  is  shown  for  the  static  pre-crush  and  the  subsequent  10" 
drop  on  the  shock  machine.  The  measured  acceleration-time  curves  for  the  10"  drop  are 
shown  in  Figure  26.  The  base  acceleration  yQ  of  the  shock  machine  table  is  compared  with 
the  acceleration  ya  transmitted  to  the  plates.  Some  protection  is  furnished  by  the  honey¬ 
comb,  but  more  could  have  been  achieved  by  using  a  honeycomb  specimen  with  a  smaller  area. 

9 .  BEAMS 

This  section  is  a  brief  summary  of  some  of  the  results  presented  in  the  doctoral 
thesis  of  S.  J.  Yim  [51], 

9.1  Test  Results  for  Beams 

Steel  beams  were  tested  in  the  arrangement  shown  in  Figure  3.  The  beams  were 
rectangular  bars  1"  x  1/8"  x  18".  At  each  end  the  beam  was  clamped  in  3"  long  clamping 
blocks.  The  beam  carried  a  central  mass  made  up  of  four  3"  x  4",  1-7/16"  thick  steel  blocks 
clamped  to  the  beam  midpoint.  The  total  weight  of  the  central  mass  and  six  bolts  and  nuts 
was  20  lb. 

The  beam  materials  were  WARPLIS  and  1018  steel.  The  chemical  composition  of 
WARPLIS  is  as  follows: 

Carbon  0.85-0.95  % 

Manganese  1.00-1.25 
Silicon  0.15-0.35 
Chromium  0.40-0.60 
Tungsten  0.40-0.60 
Vanadium  0.15-0.25 

Beams  of  1018  steel  were  obtained  in  three  conditions:  Cold-rolled,  hot-rolled,  and 
annealed  at  1200°F  for  2  hours  and  20  minutes  and  allowed  to  cool  in  the  furnace  for  17 
hours . 

Results  of  simple  tension  tests  of  the  four  beam  materials  are  shown  as  stress- 
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strain  curves  in  Figure  27.  WARPLIS  and  the  cold-rolled  1018  do  not  have  definite  yield 
points  and  exhibit  strain  hardening.  The  hot-rolled  and  annealed  1018  have  sharp  yields 
and  then  a  flat  stress-strain  curve. 

Beams  were  first  clamped  in  the  rig  of  Figure  3  and  tested  in  the  Tinius-Olsen 
testing  machine  to  obtain  quasi-static  load-deflection  and  load-strain  curves.  Typical 
experimental  static  load  versus  central  deflection  curves  are  shown  in  Figure  28  for 
WARPLIS  and  for  annealed  and  hot-rolled  1018  steels.  Then  dynamic  tests  were  done  with 
the  beams  in  the  test  rig  on  the  drop  table  shock  machine.  New  beams  were  used  in  each 
static  and  dynamic  test  so  that  the  material  experienced  first  loading  in  each  test. 

The  quantities  measured  in  the  shock  machine  tests  were:  absolute  acceleration 
of  the  central  mass,  designated  y i J  absolute  acceleration  of  the  base,  yQ;  and  strains  at 

various  stations.  Dynamic  loads  on  the  beam  were  taken  as  (mvj  +  mg).  The  relative  velo- 
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city  (yj  -  yQ)  and  the  relative  displacement  [y1  -  yQ)  of  the  beam  mass  relative  to  the 
base  were  obtained  by  numerical  integration  of  the  relative  acceleration  (vj  -  y  ). 

For  WARPLIS,  dynamic  load  versus  central  deflection  is  shown  in  Figure  30.  The 
dynamic  curve  is  higher  than  the  static  curve,  indicating  that  the  material  is  sensitive 
to  rate  of  loading  and  was  stiffer  during  the  dynamic  tests. 

Experimental  values  of  yt  and  yQ  versus  time  are  shown  in  Figure  31.  The  maximum 
acceleration  of  the  mass  is  about  twice  that  of  the  base,  showing  that  the  beam  does  not 
protect  the  mass. 

Experimental  deflection  versus  time,  strain  versus  time  and  load  versus  strain 
are  shown  in  Figures  32,  33  and  35  respectively.  These  are  compared  with  theoretical  pre¬ 
dictions,  which  are  discussed  in  the  next  section. 

9.2  Theoretical  Results  for  Beams 

The  theoretical  results  were  obtained  using  the  analytical  procedure  described 
briefly  in  Section  II. 5  and  in  detail  in  Reference  [51].  The  stress-strain  relationship 
used  is  that  of  Equation  C26) 

1  -  I  +  K(a'ast)  (26) 

Here  ojt  is  the  static  stress  existing  at  a  certain  strain  level.  In  the  analysis  vs.  e 

was  taken  as  a  bilinear  curve.  When  K  is  very  large  and  the  rate  terms  are  small,  o 

approaches  o  .  This  is  shown  in  some  typical  theoretical  dynamic  stress-strain  curves  in 
St  ,  )n2 

Figure  29  for  WARPLIS.  When  K  =  5  x  10  -  a  large  value,  the  dynamic  curve  follows 

the  bilinear  static  curve.  When  K  =  5  x  10~7  ^ys-ec,  a  small  value,  the  rate  terms  dominate 
and  the  dynamic  curve  tends  to  follow  the  elastic  curve,  -  =  !•■.  The  range  of  K  values 
between  these  extremes  produces  curves  such  as  that  of  Figure  29  for  k  =  5  x  10" 5  T-ry - , 
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the  value  which  seemed  best  for  WARPLIS.  These  curves  are  obtained  for  a  dynamic  analysis 
at  a  beam  location  one  inch  from  the  clamp  and  in  the  top  fibers,  so  the  strain  rate  varied 
considerably  during  the  loading.  The  stress-strain  curve  varies  with  strain  fate,  so  it 
will  be  different  at  each  beam  location  and  each  layer  on  the  cross-section. 

Theoretical  dynamic  load  versus  central  deflection  is  shown  for  WARPLIS  with  K  = 

.  2 

5  x  10" 5  ip .  See i n  Figure  30.  The  agreement  with  experimental  results  is  very  good.  The 
measured  base  acceleration,  yQ,  was  used  as  input  for  the  analysis.  Theoretical  accelera¬ 
tion-time,  central  deflection-time  and  strain-time  curves  are  compared  with  experimental 
in  Figures  31,  32,  and  33  respectively.  Agreement  for  accelerations  and  deflections  is 
very  good,  and  fairly  good  for  strains.  Strains  are  the  most  difficult  to  predict,  because 
they  depend  on  the  second  derivative  of  deflection.  The  location  of  strain  of  Figure  33 
was  top  and  bottom  of  the  beam,  1"  from  the  support.  The  beam  deflected  downward,  so  the 
top  fiber  is  in  tension  and  the  strain  is  positive.  The  lower  fibers  are  first  in  com¬ 
pression  (negative  strain).  The  beam  deformation  is  primarily  in  bending  during  the  first 
few  milliseconds.  However,  tension  then  begins  to  develop  due  to  the  axial  restraint  of 
the  end  clamps.  If  plane  sections  remain  plane,  then  the  axial  or  tensile  strain  at  the 
mid-fiber  would  be  the  average  of  the  strains  at  the  top  and  bottom  fibers.  It  is  seen  in 
Figure  33  that  the  average  is  zero  initially,  but  that  considerable  axial  strain  then 
develops.  The  axial  strain  is  associated  with  string-like  behavior  of  the  beam. 

In  Figure  34,  theoretical  strain-time  curves  are  shown  for  the  top  layer  and  1" 
from  the  beam  clamp  for  different  values  of  K.  At  this  station,  the  strain  is  relatively 
insensitive  to  K. 

In  Figure  35,  theoretical  and  experimental  load  versus  strain  curves  are  com¬ 
pared  at  the  top  and  bottom  fibers  1"  from  the  beam  clamp.  The  average  of  the  two  strains 
at  a  certain  load  level  would  be  the  axial  strain.  It  is  seen  that  static  and  dynamic 
curves  are  quite  close,  indicating  little  rate  of  loading  effect.  The  theoretical  dynamic 
tensile  strain  becomes  too  large,  and  bending  strain  prediction  is  less  accurate  than  axial 
strain  prediction. 

In  Figure  36,  theoretical  dynamic  load  versus  moment  at  the  support  is  shown  for 
WARPLIS  beams  with  different  values  of  K  and  different  values  of  0^  /aQ  which  is  dynamic 
yield  stress  over  static  yield  stress.  A  plastic  hinge  develops  at  the  support  but,  due 
to  strain  hardening,  the  associated  bending  moment  continues  to  increase  until  unloading 
begins  to  occur. 

In  Figure  37,  theoretical  dynamic  stress-strain  curves  are  shown  for  annealed 
1018  steel  for  different  values  of  K.  The  bilinear  static  curve  was  taken  as  elastic, 
perfectly  plastic.  A  primary  effect  of  high  rate  of  loading  on  the  annealed  steel  is  to 
increase  the  yield  stress.  If  the  material  then  remains  perfectly  plastic,  the  most  rea¬ 
sonable  representation  of  the  shape  of  the  stress  strain  curve  may  be  obtained  by  taking 
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=  1.5  0Q  and  K  =  5  x  10'1  a  very  large  value  of  K.  The  difficulty  with  this 

representation  is  that  it  does  not  accurately  predict  the  relaxation  time,  or  the  time 
delay  between  maximum  stress  and  maximum  strain. 

In  Figure  38,  the  theoretical  dynamic  moment  at  the  support  is  shown  for  annealed 
1018  steel  for  different  value  of  K  and  o^/aQ.  For  predicting  response  of  annealed  1018 
steel,  the  best  fits  occurred  for  K  =  5  x  10' 1  -j^^^with  =  1.5  oQ.  Accuracy  of  strain 
predictions  for  annealed  1018  was  poor,  due  apparently  to  the  unstable  nature  of  the  mater¬ 
ial  and  the  tendency  for  localized  flow  to  occur.  However,  accuracy  of  prediction  of  accel¬ 
erations  and  deflections  for  all  1018  beams  was  comparable  to  that  for  WARPLIS. 
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V.  SUMMARY  AND  CONCLUSIONS 


The  study  is  concerned  with  the  development  of  analytical  methods  for  design  of 
energy-absorbing  foundations  for  protection  of  equipment  in  ships  under  nuclear  weapon 
attack.  The  first  part  of  the  study  involved  looking  for  suitable  configurations  and  for 
appropriate  representation  of  material  behavior. 

There  is  a  great  amount  of  literature  dealing  with  dynamic  behavior  of  materials, 
but  more  work  needs  to  be  done,  especially  on  unloading  and  reloading  behavior  of  metals. 
However,  enough  is  known  so  that  it  is  possible  to  estimate  the  energy  absorbed  per  unit 
volume  of  material.  To  arrive  at  an  approximate  bilinear  static  stress-strain  curve,  one 
must  know  the  elastic  modulus,  the  yield  stress,  and  the  plastic  modulus.  To  extend  to 
dynamic  behavior,  the  yield  stress  and  plastic  modulus  must  be  known  as  a  function  of 
strain  rate.  The  increase  in  yield  stress  may  be  estimated  from  Equation  (1),  but  it  may 
be  necessary  to  perform  several  tests  at  different  strain  rates  to  determine  the  material 
parameters.  For  a  dynamic  stress-strain-  time  relationships.  Equations  (2)  and  (3)  may  be 
used,  again  with  the  help  of  tests  to  determine  the  parameters  for  certain  materials. 
Material  parameters,  test  techniques,  and  related  references  are  discussed  in  Section  I. 

In  addition,  the  list  of  Supplemental  References  gives  362  references  related  to  constitu¬ 
tive  equations  for  dynamic  loading  and  their  applications  to  various  structural  systems. 

With  regard  to  configurations,  axially  loaded  crushed  tubes,  frangible  tubes  and 
inverting  tubes  are  efficient  energy  absorbers  but  they  must  be  carefully  loaded  or  guided 
with  a  mandrel,  so  they  may  be  difficult  to  apply  to  ship  shock  where  the  input  may  come 
from  any  direction. 

Side-loaded  tubes  and  honeycomb  show  promise  for  use  in  ship  foundations  as  energy 
absorbers.  In  Table  2,  page  15,  it  can  be  seen  that  honeycomb  is  a  more  efficient  energy 
absorber  on  a  ft. lb. /lb.  basis.  However,  honeycomb  is  more  expensive.  Further,  side- 
loaded  tubes  could  possibly  be  straightened  in  place  after  an  impact,  and  be  re-shocked. 

Test  results  from  the  present  study  for  side-loaded  tubes  demonstrate  that  shock 
attenuation  can  be  achieved  on  the  shock  machine.  Further,  if  desired,  the  tubes  could  be 
designed  to  absorb  energy  on  two  successive  shocks.  The  dynamic  load-  deflection  curve 
rises  to  twice  the  yield  load  at  a  deformation  of  60%  of  the  tube  diameter,  as  shown  in 
Figure  12,  for  the  annealed  tubes.  Honeycomb  has, on  the  average,  a  flatter  load-deforma¬ 
tion  curve  after  a  first,  undesirable,  sharp  peak.  This  initial  peak  can  be  eliminated  by 
pre-crushing  or  dimpling  the  honeycomb.  A  further  trait  of  honeycomb,  which  may  be  unde¬ 
sirable,  is  the  slight  oscillation  of  the  load-deformation  curve  about  the  average  line  as 
the  successive  pleats  form. 

Beams  are  used  to  support  equipment  on  ships.  Also,  much  of  the  piping  on  a 
ship  deforms  laterally,  like  a  beam.  It  is  important  to  develop  dynamic  analysis  proce¬ 
dures  for  beams  to  account  for  plastic  deformation.  Some  of  the  work  completed  in  part 
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in  the  present  study  and  in  a  related  thesis  [SI]  shows  that  a  moment-curvature-time  rela¬ 
tion  ship  of  the  form  of  Equation  (7)  can  be  used  successfully  to  predict  beam  response. 

In  addition,  values  are  suggested  for  R,  n,  and  M  for  WARPLIS  and  for  1018  steel.  The 
1018  steel  was  studied  in  the  cold-rolled,  hot-rolled,  and  the  annealed  conditions.  Trans¬ 
mitted  timewise  accelerations  and  deflections  were  accurately  predicted  for  all  four  steels 
under  shock  machine  loading.  Strain  versus  time  were  predicted  with  good  accuracy  for 
WARPLIS  and  cold-rolled  1018,  but  with  poor  accuracy  for  hot-rolled  and  annealed  1018. 

For  shock  attenuation,  beams  are  not  the  best  choice,  because  the  load-deflec¬ 
tion  curve  rises  rapidly.  A  good  arrangement  would  be  a  machine  on  a  bed-plate,  protected 
by  energy- absorbing  tubes  or  honeycomb.  This  system  could,  in  turn,  be  supported  on  a 
foundation  which  is  attached  to  the  hull.  It  should  be  noted  that  the  energy-absorbing 
tubes  or  honeycomb  act  as  a  mechanical  fuse  and  limit  the  force  the  machine  applies  to  the 
foundation,  as  well  as  vice-versa.  Thus,  the  use  of  energy  absorbers  could  result  in 
lighter  foundations  and  equipment  or  it  could  lead  to  in  systems  which  could  sustain  a 
more  severe  attack. 
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Strain,  r 


Figure  1.  Various  significant  stress-strain  curves:  RB,  relaxation 
boundary;  QS,  quasi-static  stress-strain  curve;  DC,  typical  dynamic 
stress-strain  curve;  IC,  instantaneous  curve;  Dj ,  elastic  domain;  D2 , 
domain  of  possible  dynamic  stress  and  strain  states.  (from  Cristescu 
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igure  5.  Simulated  machine-snubber  configuration  using  steel 
plates  and  side-loaded  rings. 
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Figure  8.  Load  vs.  deformation  for  a  2"  long,  side- loaded  tube, 
Tinius-Olsen  machine.  (Tube  edges  cut  by  tube  cutter. 1 


Load,  lb. 


Deformation,  in. 

l'igure  9.  Load  vs.  deformation  for  four  1/2"  long,  annealed,  side- 
loaded  tubes  using  Tinius-Olsen  machine. 


Deformation,  in. 


igure  10.  Load  vs.  deformation  for  annealed,  side-loaded  tubes  using 
Tinius-Olsen  and  shock  machines. 


Load,  lb. 


Deformation,  in. 

Figure  12.  Load  vs.  deformation  for  successive  drops  1,  2,  3,  1,  li  9 
compared  with  quasistatic  curve,  for  side-loaded  tubes. 
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Time,  milliseconds 


igure  14.  Input  acceleration  y  and  transmitted  acceleration  Vj  versus 
time,  drops  1,  2,  3,  4  using  side-loaded  tubes. 


Time,  milliseconds 

Figure  15.  Acceleration  vs.  time,  vq  and  y1,  for  drop  9 
using  side- loaded  tubes  in  machine-snubber 
configuration. 
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Deformation,  in. 


figure  17.  load  rate  vs.  deformation  for  10”  drops  1  fi  2,  for  side-loaded 


tubes  . 


Load,  lb. 


Experiment 

Theory 


Deformation,  in. 

Figure  18.  Theoretical  and  experimental  load-deformation  curves 
for  10”  drops  1  ti  2,  for  side-loaded  tubes. 


Acceleration 


Figure  19.  Theoretical  and  experimental  acceleration  vs.  tine;  >'j ,  for  10' 
drop  1,  for  side-loaded  tubes. 
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Acceleration 


Load,  lb. 


800 


Figure  21.  Theoretical  and  experimental  load-deformation  curves  for 
15"  drop  1,  for  side- loaded  tubes. 


Load,  lb. 


Figure  22.  I  xperimental  load  versus  deformation  for  as -received  and  annealed 
tubes  loaded  axially. 
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Figure  23.  Hxperimental  load  versus  deformation  for  aluminum 
honeycomb . 
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30  complete  cells 


r.xperimental  accelerations,  5  ylf  versus  time  for 
aluminum  honeycomb  for  10"  drop  1  after  a  static  pre 
crush. 


WARP  LIS 


Figure  27.  Static  stress-strain  curves  for  various  steels 
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ves  for  various  steel  beams 


Stress  xlCT  ,  psi 


20 


10"  Drop 

11  msec,  programmer 


Strain  xlO3,  in. /in. 

igure  29.  Theoretical  dynamic  stress-strain  for  WARPLIS. 
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Dynamic  load  xlO 


Theoretical 
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Central  deflection,  in 


-O-  —  Experime: 


Figure  32.  Central  deflection  vs.  time  for  WARPLIS 
drop,  11  msec  programmer 
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Theoretical 


— £ —  — 4- -  Dynamic  test  (14"  drop,  11  msec 

programmer) 

— O —  —  Static  test 


-10  12  3  4 

Strain  xlO3,  in. /in. 


Figure  35.  Dynamic  load-strain  comparisons  for  WARPLIS,  14  inch  drop, 
11  msec  programmer 
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Dynamic  load  xlO 


Figure  36.  Theoretical  dynamic  load  vs.  moment  at  support  for  WARPLIS 
beam. 
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Figure  37.  Theoretical  dynamic  stress-strain  relations  for  annealed  1018  steel 


Moment  xlO'2,  in. lb. 


gure  38.  Theoretical  dynamic  load  vs.  moment  at  support  for 
annealed  1018 
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